Oscillatory behavior of a superradiating system coupled to electron reservoirs 
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We investigate a superradiating system coupled to exter- 
nal reservoirs. Under conditions where electrons tunneling 
at a rate T act like an electron pump, we predict a novel 
phenomenon in the form of oscillations with a frequency 
uj ~ y/2TT that appear in the (photon) emission intensity, 
where T is the spontaneous decay rate of a single two-level 
system. The effect, together with a strong enhancement of 
the superradiant peak, should be observable in semiconduc- 
tor quantum wells in strong magnetic fields, or in quantum 
dot arrays. 

PACS: 42.50 Fx, 42.55 Px 



Superradiance occurs in the spontaneous coherent de- 
cay of an initially excited ensemble of N two-level sys- 
tems which are interacting with a common photon field. 
The corresponding emission rate of photons is propor- 
tional to iV 2 which is abnormally large when compared 
to the incoherent decay of N independent systems. Fur- 
thermore, the emission is not exponentially in time but 
has the form of a very sudden peak on a short time 
scale ~ l/N. The phenomenon was predicted by Dicke 
in a seminal paper |Q and observed for the first time 
M in an optically pumped hydrogen fluoride gas. Since 
then, many investigations of superradiance concentrated 
on modifications through geometry effects and dephasing 
processes such as dipole-dipole interactions which had 
been neglected in the original Dicke paper ||. Since su- 
perradiance intrinsically is a many-body problem, this 
also gave the possibility to study the concept of coher- 
ence and dephasing in a many-body context Q. 

In the present Letter, we propose an extension of the 
original Dicke model that in principle opens the possibil- 
ity for an external control of superradiance. The main 
idea is to realize the superradiant 'active region' in a 
semiconductor and to couple it to external electron reser- 
voirs through tunnel barriers, thus allowing for a varying 
electron number. We predict that for reservoir conditions 
which work like an electron pump, the initial coherent su- 
perradiant peak can be strongly enhanced if the tunnel 
rate T is high enough. Furthermore, we predict a novel 
phenomenon in the form of strong oscillations of the emit- 
ted light with a frequency to that in good approximation 
is given by 
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T, there is a smooth crossover to the conventional Dicke 
peak pj in the limit T — > 0. In contrast to oscillatory su- 
perradiance in atomic systems [EL these oscillations are 
not due to reabsorption of photons, but due to tunneling 
of electrons into an active region which is described by 
many-body wave functions. The latter are characterized 
by a total pseudo spin J and a pseudo spin projection 
M for N electrons occupying upper and/or lower levels, 
including empty levels. In contrast to the original Dicke 
problem, J is no longer conserved but develops a dynam- 
ics that is driven by the tunneling process which leads to 
a coupling of J and M, whose equations of motion have 
oscillatory solutions. 

We propose a concrete realization of the model in a 
semiconductor quantum well in a strong magnetic field 
which defines degenerate energy levels for electrons and 
holes. The 'active layer' is driven in analogy to a semicon- 
ductor diode laser by injecting conduction band electrons 
and valence band holes that can radiatively recombinc. 
In contrast to the diode laser, no stimulated emission 
processes are required. 

We first describe our model Hamiltonian. A photon 
field Hp = ^2q ^QftqftQ with creation operator Oq for 
a mode Q gives rise to transitions which change an in- 
ternal degree of freedom a = (T,4) of one-particle elec- 
tronic states labeled (i,cr) with creation operator c\ a . 
The states are degenerate with respect to i with ener- 
gies £jf = — Eii = hujo/2. The electron-photon coupling 
matrix element <?q is assumed to be independent of the 
electronic quantum numbers (i, a). This condition is ful- 
filled if the momentum matrix element (icr|p|i'cr') cx 5a' 
for a ^ a' as is the case for the applications discussed 
below. Then, the Hamiltonian without reservoir coupling 
can be written as 



Hjj — TlLOoJz 



Q 
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J+ + J-)+H p , (2) 



where the operators J. 
and J z := | £\ (c\ *Ci. 



J- '■= 4,i c iA 



"]j, c i,J.) form a (pseudo) spin 



for large T > 2T, where T is the spontaneous decay rate 
of a single two-level system. For smaller tunnel rates 



algebra with angular momentum commutation relations. 
Note that so far we have not refered to the explicit form 
of the one-particle states; the Hamiltonian Eq. (Q) rep- 
resents a whole class of physical systems rather than one 
specific experimental situation. The model and the ex- 
amples below refer to a photon field; we point out, how- 
ever, that our results are valid for a coupling to an arbi- 
trary boson field (e.g. phonons), too. 
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We allow the number of electrons N in the active re- 
gion as described by Hp to vary by tunneling to and from 
electron reservoirs a = I/O ('In'and 'Out') with Hamil- 
tonians H a = J2k £ k c \ a c k,a for non-interacting electrons 
with equilibrium Fermi distributions f a . The tunneling 
of electrons is described by the usual tunnel Hamilto- 
nian H T = £ fciCTCt (*fe iCT 4,a c V + c.c.J with coefficients 
tkia t na * m general depend on the specific form of the 
tunnel barriers in real space and the one-particle wave 
functions. The total Hamiltonian is given by the sum 
H = Hq + J2 a Ha + Ht- We note that in this Letter we 
do not consider explicitly the effects of electron-electron 
and electron-hole interactions but concentrate on the ef- 
fects which result from the tunneling processes only. We 
do not expect the interactions to lead to qualitative dif- 
ferences, because in our case the carrier density is high 
and drastically varies as a function of time, which weak- 
ens and smears out the interaction effects. 

The coupling of the active region to reservoirs com- 
bines the dynamics of a (many) spin-boson problem, Hjj, 
with the transport through it. We use the master equa- 
tion approach to calculate the dynamical evolution of ob- 
servables like the emission intensity. This method, while 
only perturbative in the coupling matrix elements <?q and 
t a , has turned out to yield results that seem to be in 
good agreement with experiments both for the superra- 
diant problem |s]|5) and transport through small regions 
of interacting electrons (quantum dots) @^). We note 
that correlations between different tunneling processes 
(co-tunneling) are neglected here. 

The eigenstates of the active region are characterized 
by the total pseudo spin J and its projection M through 
J 2 \JM;{\}) = J(J + 1)| JM; {A}) and J Z \JM;{\}) = 
M\JM; {A}), where J is the total pseudo spin opera- 
tor. Here, {A} denotes all additional quantum numbers 
apart from J and M that are necessary to characterize 
the eigenstates of Hn- Radiative transitions obey the 
selection rule M — > M ± 1 with a spontaneous emission 
intensity Ijm, 

I JM = Ku rvjM, vjm-={J + M){J-M + 1). (3) 

Here, T is the spontaneous emission rate of one single 
two-level system; for radiative transition in atoms, 1/r 
is in the nano second range. 

Furthermore, the rates Tjm^j'M' for transitions be- 
tween eigenstates of Hd through electron tunneling may 
be written in the form 

r ju^rw = J2 T< * {lJM^j>M>f a (AE) 

+ Jj,M^JM(l-fa(-AE))}. (4) 

Here, T<" := (2n/h) £ fe 1%^, „,)* 5{AE - ej) is the 
tunnel rate for lead a which in fact is a matrix in (i, a) 
and depends on the energy difference AE between final 



and initial state. We neglect the energy and site depen- 
dence, furthermore the dependence on a is absorbed into 
the index a through the boundary conditions: the cou- 
pling to the reservoirs has the effect of a 'pseudo spin-up 
pump' where only <j =| electrons can tunnel in and cr =[ 
electrons tunnel out. That is, the chemical potential for 
a = I is assumed to be situated above all possible en- 
ergy differences of states \JM; {A}) with M differing by 
plus 1/2, and the chemical potential for a = O is chosen 
such that 'down' electrons can tunnel out to a — O, but 
not tunnel in. Furthermore, tunnel matrix elements for 
processes like tunneling from a — I with a =J., leading 
to M' = M — 1/2, are assumed vanishingly small; this 
condition is satisfied in the applications below. 

The coefficients 7 in Eq. (m are determined by the 
Clebsch-Gordan coefficients for adding or removing a 
single pseudo-up or down spin to the active region, 
JJM^J>M> :=\12 m=±1/ 2(J , M'\JM;j = 1/2, m)| 2 . Here, 
we adopted the same approximation that has been used 
in 0] for calculating matrix elements for transition rates 
in transport through correlated quantum dots: the spe- 
cific form of the many-particle wave function in the active 
region (i.e. the quantum numbers {A}) is neglected, and 
the matrix elements are approximated by the Clebsch- 
Gordan coefficients. 

We describe the dynamics of the active region in terms 
of probabilities p(JM) t which are the diagonal elements 
of the reduced density operator at time t in the basis of 
eigenstates \JM{X}), where all information apart from 
the quantum numbers J and M has been traced out. In 
lowest order perturbation theory, there is no interference 
between H ep and Ht, and the master equation in Markov 
and Born approximation reads S 

p (JM) t = -T{vj M p{JM) t - v JM+1 p{JM + l) t } 

— FjM->J'M'p(JM) t }. (5) 

J'M' 

The time dependence of the expectation value of the 
emission rate, Eq. (||), was obtained from the time evo- 
lution of p(JM) t by numerical solution of Eq. (||). The 
result is shown in Fig. (|l|) for an initially excited state 
with J = 6, M — 5. The initial emission maximum is 
the original 'Dicke-peak', followed by oscillations that die 
out at an intensity proportional to the tunnel rate T that 
was chosen symmetric, T = T 1 = T° . The frequency of 
the oscillations increases with y/T and follows in good 
approximation the law Eq. (|l]). In particular, the initial 
peak is strongly enhanced with increasing tunnel rate. 
This behavior is related to an initial increase of the total 
pseudo spin as can be seen from the 'phase-space' plot 
(inset of Fig. (jlj)) of the expectation values of J and M 
which both oscillate in time. 
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FIG. 1. Time evolution of the emission intensity 
Vjm = Ijm/TTujo for different transmission rates T. Inset: 
(J)t vs. (M)t for T = 64. 

These observations can be understood as follows. In 
the original Dicke problem, a fixed total number N of 
electrons with exactly one pseudo spin a on each site i is 
assumed. Here, we work in a grand-canonical ensemble 
where N varies through single electron tunneling: doubly 
occupied or empty single particle levels i become possible, 
i.e. there is no longer an exact 'half-filling'. One immedi- 
ate consequence is that a tunneling electron changes the 
quantum numbers J and M. The change J of J is pro- 
portional to M itself, j(t) = TM(t)/J(t), which follows 
considering the Clebsch-Gordan coefficients for adding a 
pseudo up spin. At the same time, M is increased by 1/2 
as is the case for out-tunneling of a pseudo down spin: 
M increases by 1 at the tunnel rate T and decreases by 
spontaneous emission at a rate Tvjm- Therefore, J and 
M obey roughly 

M(t) = ~Tv J{t)M{t) + T 
j(t)=T-M(t)/J(t), (6) 

which are governed by the two parameters T and T, the 
emission rate and the tunnel rate. After eliminating J 
in the equations of motion Eq. (^J), and approximating 
vj M ~ J 2 - M 2 , we obtain M - 2TMM + lo 2 M = 
with lu = V2TT, Eq. (§). For T > 2T, this describes 
a harmonic oscillator with frequency Eq. (0) and ampli- 
tude dependent damping. For smaller T, the oscillations 
are no longer visible and Eq. (|l|) does no longer hold. 
For T — > 0, there is a smooth crossover to the conven- 
tional Dicke peak with vanishing intensity at large times 
and without oscillations. We again point out that our 
model applies to the small sample limit of the superradi- 
ant problem: reabsorption processes of photons that may 
lead to oscillatory behavior of the intensity do not play 
any role here. 

The oscillator equations Eq. (^) in fact are the quasi- 
classical limit of the master equation Eq. (||) for J> 1. 
Neglecting fluctuations of the expectations values M(t) 



and J(i), one writes ||] the time dependent probabil- 
ity distribution as p{JM) t — 5M.M{t)&j.j{t)- The de- 
tailed form of the intensity peak and the intensity oscil- 
lations obtained in this way deviate from the exact solu- 
tion of Eq. (|5|), whereas the qualitative features remain 
unchanged which we checked numerically. In the limit of 
large J, one obtains the two equations Eq. (Q) for J and 
M. 

We now turn to the question in what physical sys- 
tems the effects described above can be observed exper- 
imentally. We note that the tunneling processes can 
be replaced with classical injection processes over po- 
tential barriers, because we have assumed that quantum 
correlation is absent between subsequent tunneling pro- 
cesses. We thus propose the system of electrons and 
holes in semiconductor quantum wells in strong mag- 
netic fields. Vertical injection of conduction-band elec- 
trons and valence-band holes into an active region acts 
like the pumping mechanism described above. In fact, 
this mechanism is exactly what is used in lasers or light 
emitting diodes with forward biased pn junctions. In our 
case, mirrors as in a laser are not necessary, in particu- 
lar stimulated emission processes must play no role. The 
strong magnetic field is necessary to have dispersionless 
single electron levels i — X, corresponding to the lowest 
Landau bands (n — 0) and guiding center X in the 
conduction and the valence bands. In this case, the in- 
terband optical matrix elements are diagonal in i. The 
correspondence with our model can be seen by mapping 
its four basic single particle states to the states of the 
electron- hole system (Fig. (|^)): the empty state becomes 
the hole (H), the pseudo-spin down electron becomes the 
empty state, the doubly occupied state becomes the elec- 
tron (E), and the pseudo-spin up electron becomes the 
state with one electron and one hole. The number of total 
electrons N in our model has its correspondence via 

N = N E + N S -N H , (7) 

where N s = $/$o is the degeneracy for a given magnetic 
flux $ (<i>o = hc/e is the flux quantum), Ne the number 
of electrons in the conduction-band, and Nh the number 
of holes in the valence-band. 

We predict that an initial optical or current excitation 
of the system leads to a superradiant peak of emitted 
light that becomes strongly enhanced if the tunneling 
rate becomes higher. Furthermore, subsequent oscilla- 
tions of the emitted light should be visible at an approx- 
imate frequency Eq. (|l|). We also expect similar oscil- 
lations to be visible as weak corrections to the injection 
current; detailed calculations will be presented elsewhere. 

As a second experimental setup, we propose an array 
of identical quantum dots, coupled to electron reservoirs 
as above. The array must have the capability to coher- 
ently radiate, where each dot has a pair of well-defined 
internal levels that allow for transitions under emission 



3 



of photons. Another possible realization could be a ge- 
ometry as in the quantum cascade laser which has been 
proposed recently as an alternative to conventional 
semiconductor diode lasers. In this case, transitions be- 
tween different electronic subbands lead to photon emis- 
sion. 



SUPERRADIANCE + RESERVOIR 



ELECTRON-HOLE SYSTEM 



FIG. 2. Correspondence of our model (superradiance with 
reservoir) with an electron-hole system. 'E' denotes an elec- 
tron in the conduction band, 'H' a hole in the valence band. 

Note that in all three cases, the photon escape time 
t = L/c (L is the linear dimension of the active region 
and c the speed of light) has to be much smaller than all 
other time scales of the problem, because in our model 
we assumed the 'small sample superradiance case' where 
reabsorption effects play no role. 

We mention that all the effects described above should 
in principle be observable not only for photons, but also 
for other bosonic fields such as phonons, or magnons. 

Finally, we point out that so far we have not addressed 
the question of phase coherence. In fact, the conventional 
superradiance is a transient process that occurs only on a 
'mesocopic' time scale with an upper boundary fll|| given 
by a phase coherence time . Inelastic processes such as 
dipole-dipole interactions jlj in general destroy the phase 
coherence between single particle states and the descrip- 
tion using the Dicke states with well-defined J and M 
becomes void. On the other hand, coherence between 
states with different J and different M is not required in 
our formalism, on the contrary this would require consid- 
eration of coherent tunneling which is beyond the scope 
of our approach. We assume that dcphasing processes 
are weak such that the time to observe the initial Dicke 
peak and some cycles of the subsequent oscillations is 
still shorter than t^. An ideal case would be a strong ini- 
tial excitation to a high initial pseudo spin J, and a high 
tunnel rate T that yields fast oscillations. Furthermore, 
strong magnetic fields in general suppress scattering rates 
although at the present state we can give no quantitative 
estimates for r^. 
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